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Abstract. We obtain expUcit bright one- and two-sohton solutions of the integrable 
case of the coherently coupled nonlinear Schrodinger equations by applying a non- 
standard form of the Hirota's direct method. We find that the system admits both 
degenerate and non-degenerate solitons in which the latter can take single hump, double 
hump, and fiat-top profiles. Our study on the collision dynamics of solitons in the 
integrable case shows that the collision among degenerate solitons and also the collision 
of non-degenerate solitons are always standard elastic collisions. But the collision of a 
degenerate soliton with a non-degenerate soliton induces switching in the latter leaving 
the former unaffected after collision, thereby showing a different mechanism from that 
of the Manakov system. 
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1. Introduction 

Recently there has been considerable interest in studying the dynamics of 
multicomponent solitons/solitary waves in view of their wide range of applications 
encompassing science and engineering [HE]. In the context of nonlinear optics, 
simultaneous propagation of multiple optical pulses or beams in nonlinear media is 
governed by a class of multicomponent nonlinear Schrodinger (NLS) type equations 
which is non-integrable in general. These multicomponent NLS equations fall into 
two categories, namely incoherently coupled NLS equations and coherently coupled 
NLS equations [1]. The integrable as well as non-integrable incoherently coupled NLS 
equations have been well studied in the literature [Il[3]. Particularly, the studies on 
integrable Manakov system [1], a two component nonlinear system with incoherent 
coupling, and also its integrable N-component generalization |1H7], have revealed the 
fact that the bright solitons of these systems exhibit interesting collision scenario which 
is not possible in their single component counterparts. This collision behaviour has been 
exploited in the construction of logic gates based on optical soliton collisions [Hill] and 
also such collisions lead to the possibility of multi-state logic [TlflO]. 

The set of coherently coupled NLS systems is another interesting class of nonlinear 
evolution equations for which much attention is yet to be paid. The term coherent 
coupling here stands for the dependence of coupling on relative phases of the interacting 
fields. A fairly general governing equation for coherently coupled orthogonally polarized 
waveguide modes in the Kerr medium (see for example. Sec. 9.4.1 in ref. [1]) is 

%,z + 5?i,TT - ^qi + {\qi? + (y\q2\^)qi + Xqlql = 0, 

iq2,z + Sq2,TT + /^g2 + {cr\qi\^ + \q2\^)q2 + Mlql = 0, (l) 

where Z and T are the propagation direction and the transverse direction, respectively, 
qi and q2 are slowly varying complex amplitudes in each polarization mode, n is 
the degree of birefringence, a is the incoherent coupling parameter and A is the 
coherent coupling parameter. Similar equations also arise in the context of short 
pulse propagation in weakly birefringent Kerr type nonlinear media [HE], where the 
co-ordinate T corresponds to retarded time. In general the system ([1]) is non-integrable. 
An integrable non-dimensional coherently coupled NLS equation closely associated with 
equation ([1]) can be written as 

iqiz - qitt - 7(kiP + 2|g2ngi - iqhl = o, 

iq2z - q2tt - 7(2|gir + 192^)^2 - 7?i?2 = 0. (2) 

The above set of equations results from equation ([T]) for the choice = (low 
birefringence limit), a = 2, A = 1, a choice which is possible in a cubic anisotropic 
nonlinear medium where the parameters A and a can be chosen separately but their ratio 
is fixed as ^ = 2 [HE] and by performing the transformations T — )■ a/t? t and Z — )■ —72;, 
where 7 > 0. Although the physical conditions for the above choice are stringent to 
obtain, we hope the exact results reported in this paper will serve as potential candidates 
in further analysis of the non-integrable coherently coupled NLS equations ([T]), which 
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have received attention recently [TTHH] . 

Motivated by the above considerations, in this paper we have obtained general 
soliton solutions of system (|2]) by applying a non-conventional form of Hirota's 
bilinearization method |15]. Another integrable equation which can also be obtained 
from equation ([1]), having a form similar to equation ([2]) but with the replacement of 
the '— ' sign appearing before the coherent coupling term by a '+' sign with 7 = 1, 
has been studied in refs. [111112] and special one- and two-soliton solutions with less 
number of parameters have been obtained by applying the Hirota's direct method. In 
fact, while obtaining the two-soliton solution by a linear superposition as reported in 
ref. [12], following the lines of ref. [11], the governing equation gets decoupled into 
two independent NLS equations and so the information regarding the coherent and 
incoherent coupling terms gets lost. This system can also be studied by applying a 
similar method as developed here and the results will be published separately. The 
main objective of this paper is to obtain an appropriate bilinear form of equation ([2]) 
resulting in more general soliton solutions as done in ref. [16] for the Sasa-Satsuma higher 
order NLS system. We also wish to investigate the soliton formation and propagation 
due to the combined effects of self phase modulation (SPM), cross phase modulation 
(XPM) and coherent coupling between the copropagating fields. Our study shows that 
there exist two distinct type of solitons, namely degenerate and non-degenerate solitons, 
where the non-degenerate solitons can have single and double hump profiles. Their 
collision behaviour is also fascinating. Particularly, the collision between degenerate 
and non-degenerate solitons shows a different kind of switching mechanism in the two 
component system ([2]) from that of the shape changing collisions occurring in the 
Manakov system [5l[7j. 

This paper is organized in the following manner. The non-standard way of obtaining 
the bilinear equations of the integrable system ([2]) by introducing an auxiliary function 
is discussed in section 2. The general one-soliton solution is obtained in section 3 and 
the degenerate and non-degenerate solitons are discussed. In section 4, the more general 
two-soliton solution reflecting the effects of coherent coupling terms during collision is 
obtained. The collisions of degenerate solitons and non-degenerate solitons are discussed 
separately in section 5 and we have also analysed the collision of a degenerate soliton with 
a non-degenerate soliton in the same section. Final section 6 is allotted for conclusion. 

2. A non-standard bilinearization method for the coherently coupled NLS 
system 

The soliton solutions of system (j2]) can be obtained by applying the Hirota's 
bilinearization method [^, which is a powerful tool for integrable nonlinear partial 
differential equations. To obtain the correct bilinear equations, resulting in more general 
soliton solutions displaying the effects of SPM, XPM, and coherent coupling, we adopt 
a non-standard method by introducing an auxiliary function, similar to the technique 
followed by Gilson et al [16] for the higher order NLS system. By performing the 
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bilinearizing transformation 

g , h 

gi = - and q2 = -, (3) 

to equation (|2]) and introducing an auxiliary function s, we obtain the following set of 
bilinear equations, 

Dig-f = -^sg\ D,h-f = ^sh*, (4a) 

D,f-f = 2^{\g\' + \h\'), sf = g'-h\ {Ah) 

where Di = iDz — D^, D2 = Df, g and h are complex functions, while / is a real 
function, * denotes the complex conjugate and the Hirota's bilinear operators and 
Dt are defined as [15] 

Note that the necessity for the introduction of an auxiliary function s{z, t) becomes 
crucial as otherwise in the absence of s in equation (4), only special cases of even one- 
soliton solution reported below will be obtained and for higher order solitons severe 
constraints on the soliton parameters will arise. The above set of equations (4) can be 
solved by introducing the following power series expansions for g, h, /, and s: 

9 = X9i + X^93 + • • • , h = xhi + x^h + • • • , (6a) 
/ = 1 + x'/2 + xV4 + . . . , s = x's2 + x's, + . . . , m 
where x is the formal power series expansion parameter. The resulting set of linear 
partial differential equations, after collecting the terms with the same powers in x, can 
be solved recursively to obtain the forms of g, h, /, and s. 



3. Bright one-soliton solutions 

In order to obtain the one-soliton solution, unlike in the Manakov case [316], here 
we restrict the power series expansion (6) as g = X9i + X^fi's? h = x^i + X^h^, 
/ = 1 + + X^fi) -s = X^'S2- After introducing this series expansion in equation 
(4) and by solving the resulting set of linear partial differential equations recursively, 
one can obtain the explicit one-soliton solution as 

^2 ~ 1 _^ ^m+Vl+Rl + g2j7i+2r?*+eii ' 

where the auxiliary function takes the form 

s={al- PDe^^'K (7c) 



Here 



e-^2fMz^, (7.) 
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Case(i): aj-^l = 

This choice — = always results in the standard "sech" profile for the bright 
soliton solution (7). It can be expressed as 



ai _Ri 



sech r]iR + 



Ri 



Aisech r]iFi + 



^1 



and q2 = ±qi corresponding to f3i = ±ai so that \qi\ 



Ri = log 



27|ai 



, and rji = r]iR + ir]ij, where rjiR 



\q2\ ■ Here Ai - 
-- kijiit + 2kijz) and rjn 



,(fei+fc*)2 

kijt + {klj — k\^z. Throughout this paper the subscripts R and / represent the real and 
imaginary parts, respectively. We call the solitons arising for the choice — = as 
degenerate solitons, owing to the fact that such solitons posses the same intensity profile 
in both the components qi and q2 and are characterized by two complex parameters 
ai and ki or four real parameters, instead of the six real parameters in the Manakov 
case [7]. Here Ai, —2ku, and are the amplitude, velocity, and central position of the 
sohton, respectively. Note that ^ is related to the polarization of the pulse/beam. The 
degenerate soliton having a single hump profile is depicted in figure 1 for the parameters 
J = 2, ki = 1 + i, and cti = /3i = 1 at t = 0. 
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Figure 1. Degenerate one-soliton a.t t = (parameters are as given in the text). 

Note that in the present case, since af — Pf = 0, the auxiliary function s vanishes, 
see equation (7c), and so from the bilinear equations (4) one can easily infer that the 
effect of coherent coupling vanishes. 
Case(ii): - /3f ^ 

The nature of soliton for the other choice af — l^f ^ 0, can be understood by 
rewriting qi and q2 in the expression (7) as 

_ 2Aj [cos(P,)cosh {rjm + + ^ sin(P,)sinh {r],R + ^)] e^''" . _ 

where A, = e(^^^^), A2 = eC-"^^) , P, = A = ^^^i^, L = 

e(^i-^) - 2, r]iR = kinit + 2kiiz), rju = kut + {kjj - A;?^)z, h = In(ai), and 



Coherently coupled bright optical solitons and their collisions 



6 



I2 = ln(/3i). Also, the quantities ^n, pn, Ri, and en are as defined in equation (7). Here 
Aj represents the amphtude of the sohton in the j-th component and for this case by the 
term amphtude we mean the peak value of the soliton profile. The speed of the soliton 
is given by 2kii and its central position is ^J^. It can be noticed that by rewriting 
the expression (9) for the choice ai = it reduces to equation (8). The general form 
presented here will be of use in the asymptotic analysis of the two-soliton and multi- 
soliton solutions. We refer to the above soliton as non-degenerate due to their distinct 
intensity profiles in the qi and q2 components. In contrast to the degenerate solitons 
these solitons can vary their profile from a single hump to a double hump through a flat- 
top profile as the parameters are varied. A double hump soliton and a flat-top soliton 
appearing in the qi and q2 components, respectively, at t = are shown in figure 2 for 
the parameters 7 = 2, fci = 1 + ai = 0.7114, and f3i = 1. Similar kind of fiat-top 




-4-20 2 4 

z 



Figure 2. A non-degenerate soliton at t = 0: (a) Double hump non-degenerate 
soliton in the qi component, (b) Flat-top non-degenerate soliton in the q2 component 
(parameters are as given in the text). 

structures have been reported in complex Ginzburg-Landau equation [2J. Note that in 
equation (9) also the standard sech type soliton occurs for a particular choice of the 
parameters, namely ai/^^ -|- al(3i = 0, for which Pi, P2, and L become zero in equation 
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(9). However, in the present case, the effect of coherent couphng does not vanish unlike 
the case of degenerate solitons. 

It can be noticed that equation ([2]) is embedded into the matrix NLS equation which 
is integrable via Inverse Scattering Transform (1ST) method [ITKIB]. Cases (i) and (ii) 
have also been reported in refs. [TTllTHj for a three component version of the equation 
considered in this paper by applying the results of the 1ST method for the matrix 
NLS equation and the corresponding solutions were referred as ferromagnetic and polar 
solitons, respectively, in the context of multicomponent spinor condensates. Here we 
have obtained similar kind of more general soliton solutions for the two component case 
itself by applying a non-standard type of Hirota's bilinearization method. 

4. Bright two-soliton solution 

The two-soliton solution of the system ([2]) can be obtained after terminating the 
power series (6) as 51 = X9i + X'gz + X'Qb + X^Qi, h = xh + X^h + X^h + x'^h, 
/ = 1 + X^f2 + xV4 + xVe + xVs, s = x^S2 + x'^Si + x^sq and again by solving the 
resultant linear partial differential equations recursively. Then the explicit form of the 
two-soliton solution can be written as 

gi = -^, J = 1,2. (lOa) 

The functions A^^^), A^^^) ^nd D in (10a) are given by the expressions 

_j_ qVi+Vi+V2+Si _|_ g»?2+'75+r?i+<52 _j_ g2j7i+2j7*+r?2+Mii _|_ g2-»;i+2'/7*+r?2+Mi2 
_|_ g2772+277*+r?i+/X2i _|_ g2»72+2»7|+r?i+/^22 _|_ g2»7i+r?f+J72+r?2+Mi 
_j_ ^2ri2+V2+Vi+Vl+fJ-2 _|_ g2r?i+2r?*+2r?2+??2+'^i _|_ g2»7i+2»72+2r?5+»?i+02 (IQ^) 
A^(2) = PiC'^'^ + 132^^^ + e^'^^'*'''!'^''" -|- e^''^"''''2+Pi2 _|_ g2r?2+»75'+P2i _|_ g2»72+»?5+P22 
_|_ QVi+Vt+V2+Pi _|_ g'72+r?2+»?i+P2 _|_ g2r?i+2»;*+r/2+i/ii _|_ g2»7i+2»7|+r?2+i^i2 
_j_ g2772 +2775" +r?i +1/21 _j_ g2»72+2r?2+»7i+i^22 _j_ ^2rii+ri^+ri2+ri2+ui 
_|_ Q2ri2+V2+m+-nt+i^2 _|_ g2r?i+2j7*+2»72+r?*+Vi _|_ g2»7i+2r?2+2r?*+»7*+?/)2 (IQc) 
£)_]__!_ ^m+Vt+Ri _j_ qVi+V2+^0 _j_ QV2+Vt+^o -)_ q'12+V2+R2 _j_ g2?7i+2?7*+eii 
_|_ g2'?i+2»?2+ei2 _j_ g2r?2+2j7*+e2i _|_ g2r;2+2r?J+e22 _|_ ^2r]i+r)l+rj*+Ti 
_j_ g2'7i+r?i+-»)2+-ri _|_ ^2'ri2+ri*+ri^+T2 _|_ ^2rf^+rii+-q2+r2 _j_ g-'?i+»7j'+»?2+r?2+-R3 
_|_ g277i+277*+r?2+»75+fli _|_ g2r?i+2j75+r?2+r?*+6»i2 _j_ g2»72+2j7*+»7i+?7j+e2i 
_|_ g2??2+2?75+r?i+»;*+e22 _|_ g2(r;i+»)*+7;2+)?|)+-R4 (10(i) 

and the auxiliary function s is determined as 

S = (a2 _ ^2)e2^i + (g,2 _ ^2^g2r,2 ^ 2(aia2 - /3l/32)e''^+''^ + e''l+''l+2r,2+An 
_|_ g'?i+^2+2^2+Ai2 _j_ g'72+f?J+2»7i+A2i _|_ gr?2+'?5+257i+A22 _j_ g2r?i+2»7*+2»72+Ai 
_|_ g2r;i+2?72+2?7|+A2 _|_ ^2r]i+'ri*+2'r)2+V2+>^3 (lOc) 
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Here 77^ = ki{t — ikiz), i = 1,2. The real and imaginary parts of rjj are given by 
VjR — kjR{t + 2kjiz) and rjji = kjit+ {k'jj — k'^jj^z, j — 1, 2. Various quantities appearing 
in equation (10) are given in the Appendix, as they are rather lengthy expressions. In 
order to understand the structure of the above two-soliton solution, we now perform 
an asymptotic analysis and analyse the nature of the soliton collisions in the present 
system. 



5. Collision of solitons 



The two-soliton solution obtained in the previous section represents the interaction of 
two solitons. It is of interest to consider the collision among non- degenerate solitons and 
degenerate solitons, and also the collision between the non-degenerate and degenerate 
solitons. For this purpose we perform the asymptotic analysis of the two-soliton solution 
(10) by considering the case where km, k2R > and ku > k2i, without loss of generality. 
The analysis is straightforward for the other choices of kjR and kji, j — 1,2. 

5.1. Collision of non- degenerate solitons (a^ ^ j — 1,2) 

The asymptotic forms of and 5*2 before collision {z — > —00) and after coUision 
{z — >■ -|-oo) can be deduced from equation (10) as follows. The quantities rjjR and rjji, 
j — 1,2, appearing in the following asymptotic expressions are defined below equation 
(10). 

1. Before Collision (z — >■ —00) 
Sohton Si {rjiR ~ 0,rj2R -00): 

ql- \^±(a\- \ ( cos(Pi) isin(Pi) \ ( coshirj-j , _ 
ft - dA Al- [ cos(P2) i sm{P2) ] I sinh(r/r«) ' ^ ' > 



where 



\ ^ _CR4±f22) 



^1- y = 2e 2 ( I , (11&) 

Di^A cosh2(r/-^) + eV'' —) - 2. (11c) 

In the above, Pi = ta^m ^ ^ t2i^^^ ^-^ ^ ^ Ra^_ Here and in the 
following the superscript denotes the soliton and the subscript denotes the component 
and - (+) sign appearing in the superscript represents the asymptotic form of the soliton 

before (after) interaction. 
Soliton 5*2 {ri2R ~ Q,riiR 00): 

- D2[ Al- [ cos(g2) i sin(g2) i I smh{rj^^) ' ' ' ^ 



where 
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2e-^ ^„-,„. , (126) 



Al~ j \ _ '2+/22 



e 2 



1^2 = 4 cosh'(r72-^) + el^2- 2 j _ 2. (12c) 

Here, Qi = ^22^, Q2 = ^22i_:^, ^ in(a2), l' = ln(^2), and 7/2}, = rj2R + ^. All the 
quantities appearing in the above asymptotic expressions (11) and (12) are defined in 
the Appendix. 
2. After Collision: 

The asymptotic expressions after collision are similar to those of before collision 
expressions with the replacement of Aj^ and rjjj^ by Aj^ and 77^, respectively for the 

sonton bj, j - wnere - (^.+fc2)(fci-fc2) ' ~ {ki+k2)(ki-ki)'^i ' ^ - -L'^^ 
r]'l^ = r]iR + and = 772^ + i^i^^iil. The quantities en and are given in 
the Appendix. One can easily check that the intensities before and after interaction 




Figure 3. Elastic collision of non-degenerate solitons (parameters are as given in the 
text). 
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are same, that is, p = l^/^P) jJ = 1)2. Also, the velocities of the two coUiding 
solitons S'l and 5*2 arc exactly the same before and after collision except for a phase 



shift which is found to be $i = = -^-In 

^ 4kiR kiR 



(fc2+fci)(fcl+fc|) 
(fcl-fc2)(fei*-fe5) 



, for the soliton and 



the sohton 5*2 experiences a phase shift $2 = —^1 ■ Thus our analysis on the non- 
degenerate sohtons arising for the general choice q;| — 7^ shows that these type of 
solitons always undergo standard elastic collision in the coherently coupled NLS system 
(2) and one such collision is depicted in figure 3 for the parameters, 7 = 3, fci = 1.5 + i, 
/c2 = 2 — i, «! = 1, /?i = 1.7, a2 = 1, and /?2 = 2. In figure 3, the double hump solitons 
undergo elastic collision in the qi component and in the q2 component the single hump 
solitons exhibit elastic collision. One can also have the double hump solitons in both 
the components, for suitable choices of parameters. 

5.2. Collision of degenerate solitons 

The degenerate solitons arise for the choice — /3| = 0, j = 1, 2. This happens when 
«! = ±/3i and 02 = ±/32. In the following, we perform the analysis for the case ai = (3i 
and a2 — ^2- For the other choices, that is, cki = ^1 and a2 — —^2 or cti = — /3i and 
Q!2 = (^2, also the collision scenario is similar to the choice discussed in this subsection. 

1. Before Collision (z — > —00) 
Soliton Si. 

q\- = ql- = Al-sech(r;^K)e^''^^ (13) 

where = ^ and 77^"^ = rjiR + . 

Soliton 5*2: 

qf- = qt = A'-sech{ri^j,)e'^^', (14) 

where = 2^ and 77^^ = rj2R + 

2. After Collision (z ^ -\-oo) 
Soliton Si. 

gJ+ = g^+ = Ai+sech(r7+,)e^''-, (15) 

where — ^ 77^ = r]iR + 
Soliton 52: 



where A^^ = ^ and = r)2R + ^ 



gr = ?2^ = >l'''sech(r/+Je^''-, (16) 

and ntR = r]2R^^^. 
All the quantities appearing in the above expressions (13-16) can be obtained from 
the corresponding quantities defined in the Appendix with the substitution /3j = a^, 
J = 1,2, and the real and imaginary parts of ?7j-s are defined below equation (10). 
Prom the above expressions, one can show that the amplitudes A^-% before and after the 
interaction are related through the expressions A}'^ 



{k*-k*)(kl+k2) 



A^- and = 

which shows that the intensities before and after interactions are 



{k*~-k^){ki+k^) 
(ki-k2)ikl+k2) 
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the same, that is = \ A'~\'^,j — 1,2. Also the sohton 5"! undergoes a phase shift 

$1 = ^^^kiR^^ ' whereas the sohton 5*2 experiences a phase shift $2 = —^1 during 
cohision. Thus the degenerate sohtons always undergo standard elastic collision as that 
of the NLS solitons. 

5.3. Collision between degenerate and nan- degenerate solitons 

The collision of a degenerate soliton (q;| — /3j — 0) with a non-degenerate soliton 
{a'j—l3j 7^ 0) exhibits very interesting collision properties. Here we consider the collision 
of a non-degenerate soliton Si (ai 7^ with a degenerate soliton S2 {0^2 = 1^2) ■ Note 
that the analysis can also be performed for the other possible choices like a2 = —(^2, 
but here also one can infer the same kind of collision scenario as for the present choice 
a2 — /32- The asymptotic forms of the sohtons Si and S2 are presented below. 

1. Before Collision 
Soliton 5*1: 

cos(Pf ) i sin(Ff ) \ / cosh(?7f^) \ ^^^^^ 
cos(gj") i sin(Q^) J y smh{r)^j^) J 

■^2+Ml-9ll--«2 \ 
6 2 / 

+ ^11^, and L^- = eV"^ 2 ; _2. Note that the expressions 
for various quantities appearing in equation (17) and in the equations (18-20) given below 
can be obtained from the corresponding quantities defined in the Appendix by putting 
(^2 = a2. 
Soliton 5'2: 

qf- = qj- = ^^-sech(r;-^)e"'-, (18) 

where A^' = ^ and ry^ = r]2R + 

2. After Collision 





Soliton ^i: 




C0S(P+) 2sin(P+) \ / COSh(7]+j) \ j^^^ 

cos(Q+) I sin(g+) i I sinh(r7+,) ' ' ^ ^ 



= 4cosh2(77+,) + /+ = In(ai), l^j = ln(/3i), 
ViR = ViR + = ~ ~ 2. The amplitudes A]^ and Al'^ are given by the 

relations A\'^ = Ti A\~ and A^'^ = T2 Al' . Here the transition amplitude 

' \ \[{ki-k2f + {ki + klf]ai + {k2-kl-2ki){k2 + kl)fii\^ 
and the expression for T2 can be obtained by replacing ai •H- j3i and a\ •H- in the 
expression for Ti. 
Soliton S2'- 

ql-^ = ql+ = A2+sech(772+K)e^''- , (21) 
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where A^^ = and = ri2R + The real and imaginary parts of 

T]j-s appearing in the above expressions are aheady defined below equation (10). 

From the asymptotic analysis we observe that the amplitudes and hence the 
intensities before and after collision are not same for the non-degenerate soliton ^i, 
while it is so for 82- It should be noticed that the arguments of the circular functions are 
also different before and after collision and also 7^ L^+. Additionally, there occurs 
a phase shift, $1 = , for soliton Si. Thus, during its collision with soliton S2, 

the soliton 5*1 experiences an intensity switching among its two components resulting 
in a redistribution of the amplitude and phase. But the amplitude of the other soliton 
5*2 remains unaltered as = and hence it undergoes the standard elastic 

collision only along with a phase shift $2 = ^""^^"^"'^^ ■ However, ^2 induces the collision 
with shape changes (intensity redistribution) in soliton Si during collision. This collision 




Figure 4. Shape changing collision of a non-degenerate soliton with a degenerate 
soliton S2 (parameters are as given in the text). 

scenario is quite different from the shape changing collision occurring in the Manakov 
system [5H7], where there is an intensity redistribution among the solitons in both the 
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components but in the present system it happens only among the two components of 
the non-degenerate sohton 5*1. Note that though the total energy of both the solitons 
is conserved independently due to the conservation law, + |q'2p)c^^=constant, 

the energy of the soliton in the individual modes that is \qi\'^dt and \q2\'^dt 
are not conserved independently as the soliton 5*1 only experiences intensity switching 
in both the components. It could be an interesting future study to check whether \Tj\, 
j — 1,2, can be unimodular, if so, for what choices of a-s and ^-s this will happen. 

For illustrative purpose, the above collision scenario is shown in figure 4 for the 
parameters j — 2, ki — 2.3 + i, k2 — 2.5 — i, ai — 0.75, Pi — 1.9, and a2 — — ^ + i- 
The figure shows that in the qi component the single hump soliton changes its profile 
to a double hump soliton and also experiences significant suppression in its intensity 
whereas the soliton undergoes elastic collision. The reverse scenario takes place for 
the soliton Si in the 52 component and here also the soliton 5*2 remains unaltered after 
collision. 

In the collision of non-degenerate solitons alone the coherent coupling modifies 
uniformly both the solitons before and after collision, thereby resulting in an elastic 
collision. But in the present case the effect of coherent coupling is switched off in the 
degenerate sohton 5*2 (since q;| = ^|), however the coupling still persists in the non- 
degenerate soliton Si. Hence, along with the XPM term the coherent coupling influences 
the non-degenerate soliton Si resulting in an intensity switching during collision. From 
a mathematical point of view, one finds that the asymptotically dominant terms of the 
non-degenerate soliton collision case become insignificant and the less dominant terms 
in that case become significant in the two-soliton solution expression corresponding 
to the collision of a degenerate soliton with a non-degenerate one. This yields different 
asymptotic expressions for these two collision processes as seen in the present subsection 
and in section 5.1, which ultimately makes their collision scenario completely different. 

6. Conclusion 

Explicit forms of one- and two-soliton solutions of the coherently coupled NLS equations 
have been obtained using a non-standard type of Hirota's bilinearization method. 
Analysing the nature of the bright one-sohton solution we have reported degenerate 
solitons (solitons possessing same intensity in the qi and 52 components) and non- 
degenerate solitons (solitons with different intensities in the qi and q2 components). 
Particularly, for non-degenerate solitons the density profile can vary from single hump 
to double hump profile including fiat-top solitons. Our analysis on the collision dynamics 
revealed the fact that separate collisions among degenerate solitons alone or among non- 
degenerate solitons alone arc clastic. On the other hand, collision of a degenerate soliton 
with a non-dcgcncratc soliton exhibits nontrivial behaviour resulting in an intensity 
switching of the non-degenerate soliton spread up in the two components leaving the 
other soliton unaltered. This property will have immediate applications in soliton 
collision based computing. Apart from the switching, we have also observed that this 
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collision transforms the soliton profile from single hump to double hump including flat- 
top profile or vice versa. We expect that this property can find application in pulse 
shaping in the context of nonlinear optics. 

The above analysis can be extended to the study of three and higher order soliton 
solutions. The details of multi-soliton collisions and the multicomponent cases will be 
pubhshed separately. 
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Appendix 

The various quantities occurring in equation (10) and in section 5 have the following 
forms: 

g,, _ a*{aj - /3f)-f _ {a*{aia2 - /3i/32)j + {h - k2){aiK2j - a2K-ij)) 



K% + k*nks-, + k^k, + k*nk,., + k*y I ^ 
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(k, + k;){ks., + k*) 

{k, + k*y {ki + klY {k2 + kiy 

^\k,-k2)\kl-kina\-Pl){al-^l) ^ . . .2 .*2^ 

2(A;, + fc*)2(fc, + A;*_^.)^ ' 2{k, + k*^\k]^k^.,Y ' 

^ (/Cl - k2fl\a\ - f^i) ^^^^^ ^ ^,^2 ^ _ ^ 

-(A:* - kl){k2 + fc2*)«2/92/5r + {k2 + - kl)a*Ml - {k^ + A;t)(fc2 + k^a^^l^l) , 

{h - k2)WH - ^ ^*)2 ^ _ ^ k,)]a,alal 

-{ki - k;){k^ + k;)a;^i^; + {k, + ki){ki - k;)aii3ip; - {h + k*^){ki + kDa^pip;) , 

-{k,-k2f^\al-Pl) ^^^^^ ^ ^ _ ^ ^,)]^,^*^* 

-(/c* - kl){k2 + A;2*);52*«2«t + {k2 + - k*)l3la2a* - {k2 + A;*)(fc2 + kl)^2alal) , 

-(^1 - k2){ki + A;*)/?*^!^! + + kl){kl - kl)/3laia* - {ki + kl){ki + k*)/3iala2) , 
{kl + fc2)^(fci + k*^^ - {ki + + k2){ki + A;2*)(A;2 + fep + (fci + fci*)^(A;2 + k*^)^ 

(ki + fci*)2(A:i* + A:2)2(fci + A:2*)2(fc2 + k*^^ 
{ki - k2){kl - A;*)(«2a2/?i/3i* + «i«t/32/32*: 



eM2 ^ 



e«3 ^ ^2 



(q;iQ;*q;2Q;2 + P1P1P2P2) + 



(k, + A;t)2(A:t + k2)ik, + A:2*)(A;2 + A:!)^ 
ala2Pif32 + aia2/3lP2) \ f {ki - k2){kl - k2){aia2/32Pl + a2alPi(32 



Q-.J — 



{ki + kl){kl + A;2)(A;i + A;2*)(A;2 + kl) J \ (ki + kl){kl + k2Y{ki + klf{k2 + k^] 
{h - k2f{kl - kl)\a^ - pl){af - /3f )(«3-i«3-,- + ^^-i^l-^W 



2D{ki + k*Y 



{ks-. + A;*)2(A:, + k*) ' 4(A:,- + A;*)4(A;; + k^.j)'^ 

where 

D = 2(A;i + A;*)2(A;t + A;2)'(A;i + kl)\k2 + A;;)^ 
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